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Abstract 

We demonstrate that p-adic analysis is a natural basis for the construction 
of a wide variety of the ultrametric diffusion models constrained by hierar- 
chical energy landscapes. A general analytical description in terms of p-adic 
analysis is given for a class of models. Two exactly solvable examples, i.e. 
the ultrametric diffusion constraned by the linear energy landscape and the 
ultrametric diffusion with reaction sink, are considered. We show that such 
models can be applied to both the relaxation in complex systems and the rate 
processes coupled to rearrangenment of the complex surrounding. 

1 Introduction 

The concept of hierarchical energy landscape attracts a lot of interest in connec- 
tion with relaxation phenomena in complex systems, in particular, glasses, clusters, 
and proteins [||]-|l2l. This concept can be outlined as follows. A complex system 



is assumed to have a large number of metastable configurations which realize local 
minima on the potential energy surface. The local minima are clustered in hierar- 
chically nested basins of minima, namely, each large basin consists of smaller basins, 
each of these consisting of even smaller ones, and so on. To be more definite, the 
hierarchy of basins possesses ultrametric geometry. Finally, the basins of minima 
are separated from one another by a hierarchically arranged set of barriers, i.e., 
high barriers separate large basins, and smaller basins within each larger one are 
separated by lower barriers. 

Transitions between the basins are determined by rearrangements of the system 
configuration for different time scales. 

Thus, two key points of the concept of hierarchical energy landscape should 
be marked. First, the configurational space of the system is approximated by an 
ultrametric space; second, the configurational rearrangements of the system are 
described by stochastic motion in the ultrametric configurational space. 

This concept has been implemented in a number of toy models referred to as 
random walks in ultrametric space, diffusion in ultrametric space, and ultra-diffusion 
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T3|| -[0. A recent term for such processes is "basin-to-basin kinetics" 0. In fact, 
all these models deal with a certain type of stochastic motion, which in this paper 
we call diffusion in ultrametric space, or shortly, ultrametric diffusion. 

The toy models have definitely demonstrated that ultrametricity actually reflects 
the characteristic features of relaxation in complex systems. It is commonly recog- 
nized that relevant analytical tools should be developed for applying the concept of 
hierarchical energy landscape to the description of relaxation phenomena in complex 
systems ^ [H) EH- In this paper, we show that p-adic analysis is an adequate 
analytical tool for studying these problems. 

The field of p-adic numbers is the most important example of ultrametric spaces 



(see, for instance, [T^])- An introduction to p-adic analysis can be found in |T8[. The 
p-adic mathematical physics attracts a great deal of interest in quantum mechanics, 
string theory, quantum gravity [|TH| - spin-glass theory |^, and theoretical 



biology [0. In works [^2] , |23| the p-adic parameterization of the replica Parisi 
matrix was developed. In papers ^ the p-adic Fourier transformation (called 
there the replica Fourier transformation) was applied in the replica method. In 
the work it was shown that the p-adic change of variables maps the basis of 
eigenvectors of the Vladimirov operator of p-adic fractional derivation onto the basis 
of wavelets in L^(R). 

We demonstrate here that p-adic analysis is a natural basis for the construction 
of a wide variety of sufficiently complex models, which can be efficiently used in ap- 
plications to both the relaxation in complex systems and the rate processes coupled 
to the relaxation of complex environment. 

This paper has the following structure. In section 2, a general analytical descrip- 
tion in terms of p-adic analysis is given for a class of ultrametric diffusion models. 



the so-called pure models (see IjT^l)- Since p-adic analysis seems to be a fairly new 
analytical tool, we thought it appropriate to describe the technique of solving the 
corresponding p-adic master equation. Specific examples are considered in the next 
two sections. In section 3, we give a complete description of the ultrametric diffu- 
sion constrained by the linear energy landscape. In section 4, we consider a special 
case when the rate process coupled to the configuration rearrangements can be de- 
scribed; and we also introduce and study a model of ultrametric diffusion with a 
reaction sink. We show that such models can be applied, in particular, to the ligand 
rebinding kinetics in heme proteins. 



2 Ultrametric diffusion 

In this Section, we describe a class of pure ultrametric diffusion models corresponding 
to regular hierarchical landscapes with degeneracy. First, we review a model of 



diffusion on an ultrametric lattice |[T3[ (this model was considered in with the 
help of p-adic analysis). Then, passing to a continuous description, we obtain the 
p-adic master equation. The last part of this section contains a brief description of 
the standard analytical technique for the investigation of such p-adic equations. 
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Ultrametric lattice and the master equation 



Consider a system with the state space Bn consisting of the points i = 1, . . . , 
where p is a prime number. These points can be regarded as lattice sites. On 
this lattice, a regular hierarchical energy landscape can be constructed as follows 
(see figure 1). Let us divide the space into p mutually disjoint subsets (basins) 
Bn-i{cli), ai = 1, . . . ,p, where each BN_i{ai) consists of p^""*^ points: U Biy_i{ai) = 

0.1 

Bn ■ 

Let us introduce activation barriers. The activation barrier for the transitions 
between the basins Bn-i{cli) with different values of ai is taken equal to Hn- Thus, 
the probability of transition between any two sites i G BN~i{ai), j G BN-i{a'i) of 
different basins Bn-i is equal to pn- In the Arrhenius case, we have pn = e~~i^ . 
Further, let us divide each basin BN-i{ai) into p smaller basins B n -2{.cii0'2) i 0,2 = 
1, . . . ,p, each having sites: lJa2 BN-2{.ciia2) = BN-i{ai)- 

Assume that the activation barriers within each larger basin BN-i{cii) are equal 
to Hn-1- Moreover, assume that the probability of transition between any two 
sites i G i?Ar_2(aia2), j G BN-2{cii0'2) of different basins Bn-2 is equal to Pn-i for 
Oi = a[, and is equal to pat for ai ^ a[. In other words, the probability of transition 
between sites of different basins Bn-2 depends on the hierarchical level, — 1 or 
A^, at which these basins merge into a single super basin. We proceed like this until 
we reach the first level with each basin Bi{aia2 . . . otv) having only one site. As a 
result, the probability of transition between any two sites is specified by a single 
value p^ among pi > . . . > > . . . > pN, according to the level 7 = 1, 2, . . . , A^ on 
which these two sites happen to be in a single basin 5^. 

Denote by fi{t) the probability of finding the system in the site i at time t. Using 
the equivalence between the transition probability matrix and the Parisi matrix 
p^, we can write the master equation for the evolution of probabilities in the 



form 22 



d 
dt 



f(t) = (Q-AoI)f(t), 



where f(t) = {/i(t),., 
matrix of Parisi type 



is the state vector, Q is the transition p x p 



[, I is the unity matrix, and Aq = J2'y=iP^~ P-y is the 
eigenvalue of the matrix Q corresponding to the eigenvector with equal components. 
For instance, the transition matrix Q for p = 2 has the form 
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This form of the transition probabihty matrix Q is a direct consequence of the 
hierarchical picture of basin-to-basin transitions. 

The master equation (0) coincides with the equation of ultrametric diffusion 
considered in ITS 



Continuum limit 



The continuum limit for the above hierarchical lattice model is discussed in 22 and 



| 23[| . The results obtained in ||22| can be summarized as follows. We introduce the 



following one-to-one mapping of integers i enumerating the lattice sites onto a set 
of real numbers x: 

n N 

3=1 i=i 

The value x*-*^ is regarded as a coordinate of the site i. Then, in the contin- 
uum limit, the set of coordinates x is extended to the field of p-adic numbers Qp, 
the matrix (Q — AqI) turns into an integral operator, and the master equation (|l]) 
transforms into the following equation: 

^f{x, t)= I {f{y, t) - fix, t))p{\x - y\p)dM . (3) 
at JQp 

Here, the operator on the right-hand side is called the ultrametric diffusion operator, 
the function p{\x — y\p): Qp x Qp H (R is the field of real numbers) is the 
kernel of the ultrametric diffusion operator, the function f{x): Qp ^ H. is the 
probability density distribution, and dfi{y) is the Haar measure on the field of p- 
adic numbers Qp. 

Expression specifies the general form of ultrametric diffusion operators. A spe- 
cial case of the operator on the right-hand side of should be mentioned, namely, 
the Vladimirov operator 

D:f{x,t) = -^ f ^^^^^^dM, (4) 
Tp{-a) Jqp \x - y\j+°' 

1 - 

where a > and Tp{a) = — is the p-adic gamma function. This operator 

is an analog of the differentiation operator in p-adic analysis. For this reason, the 
equation 

^^f{x,t) = -D:f{x,t) (5) 

is interpreted in the p-adic mathematical physics as the equation of Brownian motion 
on a p-adic line. More details concerning pseudo-differential operators and equation 
can be found in the monograph |jl8[. It is interesting to note that an equation 
of type @ is the master equation for the model in which the height of the barriers 
separating the basins of states has linear growth with respect to the hierarchy level 
number [p2 |. 
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Solution Technique 

Next, we give a brief review of the standard technique for solving equations of type 
(^. We investigate the Cauchy problem for equation (^) with the initial condition 
/(x, 0) = S{x), where 6{x) is the p-adic delta function [l^. Let us find a fundamental 
solution of equation (|). Applying the p-adic Fourier transformation to (|]), we 
obtain the following equation: 

^i%^ = -pmp)fiu), (6) 

where 

P(l^lp) = / Pi\x - y\p)il - Xi^x))dfi{x) , (7) 



and f{C,,t) is the Fourier transform of the distribution f{x,t). Hence, it is easy to 
find the Fourier transform 

/(e,t)=exp(-p(|eWt). (8) 

Applying the inverse Fourier transformation to we obtain the fundamental so- 
lution 

f{x,t)=[ exp(-p(|^|p)t)x(-ea:)rfMO , (9) 



where x{^^) is the additive character of the field of p-adic numbers. 

Equation (^) is defined on the whole field of p-adic numbers Qp. Therefore, 
the solution (|) describes the diffusion process in an unbounded ultrametric space. 
However, for some applications, it is important to study the problem of ultrametric 
diffusion in a bounded region of the state space. In this case, the following procedure 
can be used. 

Consider diffusion in the p-adic disc Br = {x : \x\p < p' } of radius p"^ with center 
at zero. In this case, the master equation reads 



df{x,t) 



{f{y,t)-fix,t))p{\x-y\p)dfx{y). (10) 



dt 

Changing the variables, x = p~'^z, we transform this equation to 

= l^i^', t) - f{z, t))p{\z - z%)dp{z') . 
This equation can be investigated by means of p-adic Fourier series 

<Piz,t) = J2Mt)x{-kz) , 

fee/ 

where the characters {x{kz)} form an orthonormal basis in L^{Zp) labeled by 
k = 0, k = p^'^ [kQ + kip + . . . + k^^ip^""^^^ 
7 = 1,2,...; fco = l,2,...,p- 1; fc^- = 0, 1, . . . ,p - 1; j = 1, 2, . . . , 7 - 1 
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3 Diffusion on linear hierarchical landscape 



In this section, we examine ultrametric diffusion in the case of a hnear landscape 
with the height of activation barriers having hnear growth with respect to the 
number 7 of the hierarchical level. Since \x — y\p = p"' , for the linear landscape we 
have 

H, = H,\Yi{\x-y\p) , (11) 

where Hq is a scale parameter. Suppose that the probability of transition between 
the states separated by the activation barrier is defined by 

= WoP"'^exp (^-^^ , (12) 

where Wq is a pre-exponential factor, T is the temperature, k is the Boltzmann 
constant. As noted in the previous section, under these assumptions, the kernel 
of the ultrametric diffusion operator coincides with the kernel of the Vladimirov 
operator (see also |2^): 



p{\x-y' ' ^° 



p 



\x - y\ 

where a = H^/kT. The Fourier transform of this kernel is 



p(0 = -^orp(-«)|e|^. (13) 

Thus, in the case of a linear hierarchical landscape (0) (under the assumption 
(HI))) the master equation of ultrametric diffusion coincides with the equation of 
Brownian motion on the p-adic line [|18| (see relations (^) and (^ in section 2): 

_ nj,Tp{-a)D:f{x, t) = . (14) 

Let us examine the solutions of this equation. As the the initial value we take 
/(x, 0) which is constant inside a bounded set in Qp and vanishes outside. For 
instance, we can take as the initial function the indicator of Zp-. 

f{x,o) = n{\x\p) = l^l'^ (15) 

Note that the p-adic Fourier transform of this indicator function is the same func- 
tion. The diagram of the model with this initial condition is given in figure 2. Using 
the fundamental solution (^, one can construct the solution of equation (|l^ with 
the initial condition (|T3|): 

fix,t)= I x(-ea:)f^(|eWexp{r,(-a)|e|>ot}ciMO- (16) 
Calculating the integral, we get 



/(x, t) = (1 - p-^) \x\p^ p~'^ I h^p 1 h^^pi-^)^ot \ - 

7=0 




\x 
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-|x|-iQ(^)exp|^r,(-aK4 • (17) 



OC If) / I X 



p 



The solution (0) describes the diffusion process in the unbounded ultrametric 
space. As an illustration, the time-dependence of population densities for different 
states is shown in figure 3. It is easy to see that at fairly low temperatures (with a 
much larger than 1), the decay of the population density of the states has a discrete 
character due to the hierarchy of the characteristic times of transition through the 
respective activation barriers. As the temperature increases, the discrete nature of 
decay becomes less and less noticeable and at sufficiently high temperatures (with 
a of the order of or less than 1) the population density of the states is governed by 
the power law. 

This conclusion can also be obtained by analytical means. Now, let us estimate 
the solution ( p!7D for \x\p > 1. Using the Abel transformations 



f{x,t) = \x\p'^J2p 



oo 

7 

P 

7=0 



exp <^-j-^Tp{-a)wot | - exp | ^^^^ Tp{-a)wot ^ 



\x\ 



+ 00 „^ J 

p ' Ep'' ^ hp {p-"1x|;"r,(-a)^ot} 



dz 



and the estimate 



we get 



pm rm rm 

p-"" / g{z)dz < / p-'g{z)dz < / g{z)dz 

Jm—l Jm—1 Jm—1 



"""^^''^ + -' -rp(-«)^ot (j^Y) < fix, t) < 



p(-rp(-«))V" ' V ' a' ' " \\x\, 

where 7(0, h) is the incomplete gamma function. If 

Wot > 



P 



\x\ 



\x\ 



p J -rp(-a)' 

then 



.(i + i-r,(-.w(j;-) )«r(i + i 

where r(a) is the gamma function. Thus, for any state x, there is a characteristic 
time 

/ I I \ a -1 
/ \ I \-^\p \ ^0 



\ p J -rp(-a) 

such that the population density f{x, t), for t ^ r{x, a), can be estimated by power 
functions: 

< H..t) < ^^(jJ-|^(u,or"" (18) 
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Let us calculate the average of the distribution (|T7|) . Note that the moment 
integrals 

{H)= I \xfJ{x,t)d^^{x) (19) 

are divergent for all /5 > a, and therefore, neither the mean p-adic distance nor 
the mean square of the p-adic distance contain any information on the distribution 
f{x,t) if, for instance, a < 1. To describe the ultrametric diffusion process in terms 
of mean values, some meaningful characteristics should be introduced. The mean 
value of the activation barrier to be overcome by the system by the instant t can be 
taken as one of such characteristics. In the case of a linear hierarchical landscape 
(see (pUP), the mean value of the activation barrier is 

{H{t)) = hJ In \xU{x, t)dfiix). (20) 

Note that the integral on the right side of ( pil| ) should be understood as the limit 

/ \n\x\pf(x,t)d^(x)= lim / In \x\pf (x,t)d^(x) . 

Let 

L^''\t)= hi\xU{x,t)d^i{x) . (21) 

J\x\p<p-^ 

By (|16D we get 

L^^){t)=f'{n{l-p-^)-p-^)\np p^exp{rp(-a)p°^u;oO- 

7=— oo 

+ 00 

- Inp exp{rp(-a)p°'^Wot} . 

7=— n+l 

Therefore, the quantities L^'^\t) satisfy the recurrent relations 

-n+l 

L^^)(p^t) = L("-^)(t) +p"-^(l -j9-^)lnj9 Y exp{Tp{^a)p''^Wot} . 

7=— oo 

Passing to the limit as n ^ +c>o, we obtain the following functional equation: 

L(p°t) =L{t)+\np 

with the solution 

L(t) = a^Mnt . 

Therefore, 

{H{t)) = kTH^Hnt. (22) 

Thus, the linear hierarchical landscape has a remarkable feature, namely, the 
mean value of the activation barrier to be overcome by the instant t has logarithmic 
growth with respect to time. 
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Summing up the results of this section, we emphasize three important features of 
ultrametric diffusion. The first and foremost is that the evolution of the population 
density in any state has a discrete nature and is determined by the hierarchy of 
the relaxation time scales. This feature stems from the concept of the hierarchical 
energy landscape and is typical for ultrametric diffusion. The discrete pattern of 
ultrametric diffusion can be observed at sufficiently low temperatures (with the 
dependence of a on T considered above). Next, in the case of a linear hierarchical 
landscape, the long-time behavior of the population density of states is described by 
the power law. Finally, the mean value of activation barriers overcome by a point 
by the time t is logarithmically increasing with the time. 

For glasses and proteins, relaxation processes with similar properties are dis- 
cussed in[|,i, i-i, [11-0. 



4 Rate processes coupled to configurational rear- 
rangements 

In this section, we apply the p-adic analysis to describe the kinetics of transforma- 
tions coupled to configurational rearrangements of the surrounding. We consider 
the well-known example of such processes, namely, the ligand rebinding kinetics of 
the myoglobin (Mb) [|, ^ |^, In this case, the object of investigation is the 

rebinding process governed by conformational rearrangements of the protein itself. 
These investigations were the origin of the concept of hierarchical energy landscape 
in the conformational dynamics of proteins. 

The experiment is as follows. The myoglobin is transformed by the laser pulse 
from the ground state [Mb-CO] into the excited unbound state [Mb*-]. Then the 
protein macromolecule relaxes into the functionally active state [Mbi-]. From this 
state, the myoglobin can go back to the ground state [Mb-CO] owing to the binding 
of the CO. 

This process can be represented as the following sequence 

Mb-CO [Mb* ^ . . . ^ Mbi] ^ Mb-CO 

The rebinding kinetics is described by the survival probability function S{t), i.e., 
the probability of the fact that myoglobin is still unbound by the time t after being 
exposed to the laser pulse. Thus, the quantity observed in these experiments is the 
total population of the conformational states [Mb* —>■...—!> Mbi] through which 
the protein macromolecule relaxes from the excited state [Mb*-] to the functionally 
active state [Mbi-]. 

Similar-type processes correlate with the models of ultrametric diffusion de- 
scribed by the master equation in the following form: 



df{x,t) 
dt 



{f{y,t)-f{x,t))p{\x~y\p)dfi{y)+ K{x,y)f{y,t)dfi{y). (23) 
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The last integral on the right-hand side of (^) can be thought of as a reaction sink 
distributed over a bounded region in the conformational space. In this case, the 
survival probability function is 

S{t) = [ f{x,t)di^{x) . (24) 



The procedure of solving equations of type ( pSD is as follows. Let us take 
f{x,0) = fo{x) in the initial condition. Consecutively applying the p-adic Fourier 
transformation in x and the Laplace transformation in the real variable t, from (^4]), 
we get 

where the symbol "marks the Laplace transform with respect to the variable s, and 
p(l^lp) is given by ([?]). Next, we note that equation ( P3| ) is a heterogeneous Fredholm 
equation of the second kind. 

Consecutively applying the inverse Laplace transformation and the inverse p-adic 

Fourier transformation to /(^,s), we obtain the solution f{x,t) and the function 
Sit). 

Let us give an example of a model of this type (see figure 4) by making the 
following assumptions: 

1. Conformational rearrangements are restricted to a bounded region of the con- 
formational space of the system and this region is the p-adic disc Br = {x : 
\x\p <p'',r:^ 1}; 

2. The hierarchical landscape is linear, i.e., conformational rearrangements are 
described by the ultrametric diffusion operator of type (H) defined on the disc 

Br, 

3. In the region Br, there is a certain set of conformational states in which the 
system can undergo irreversible transformations. This set of states we describe 
by the unit disc Zp. In other words, there is a reaction-type sink at each point 
of the disc Zp] 

4. The rate of the reaction sink at each point of the disc Zp is proportional to 
the mean value of the population density on that disc; 

5. The initial distribution f{x, 0) is constant on the the p-adic layer 

R^s = {x : p^^^ < \x\p <p'',0<6<a<r} 
and vanishes outside the layer. 
The master equation for the model with the above assumptions is 

^^=f ^^]:'l~JiV^ My)-m^\p) I ^{\y\r)f{y.t)dKy) (26) 
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and the initial consitions have the form 

/(x, 0) = AGl{x) = A [nilxlpp-") - n{\x\pp~^)] , A-^ =p'' ~p\ 0<6 < a, 

where A is the rate parameter of the reaction sink. Let us use the technique of 
Section 2 to solve equation (^61). First, we pass from the disk Br to the disk Zp by 
letting X = p~^z. Then the Cauchy problem takes the form 



dt 



p 



(t){z',t)-(t){z,t) 

ry \y yf\CL-\-l 



~ Xp'n{\z\pp') / ni\z'\pp'-)(i)iz',t)dfiiz') 

(jy{z,0) = AGrr{z), 

Using the p-adic Fourier series for Zp (see Section 2) and reasoning as above, we 

obtain a system of equations for the Laplace transforms of the coefficients of 

the p-adic Fourier series: 



k = : scpQ^s) = p — \p 



geA{o} 



A; ^ : s0o(s) = AGr_:{k)^p"'\g^ - \k\ 



-9.{\k\pp-')\p- 



qG/\{0} 



where 



9c 



1-p-^ 

1 _ 



, Grz{^)=p'^-^^{\k\pp'^-^)-p'^^^{\k\pp'-^) 



J>—T\ 



It is easy to find the solution to this system and write the equations for 0o('^) 
and If we note that in this case the following equations hold: 

S{t)=f f (P{z,t)dfi{z)=p^Mt) , 



S{s) =p''(j)Q{s) , 

we can find the Laplace transform of the survival probability function. 



S{s) 



s + \p ^ 
XA 



l + X'p 



2„-2r 



Jr(s) 



S + Xp-'-{l + Sjr{s)) 
p^-^Jr-As)-p'-^Jr-5is: 



S + Xp-'^{1 + Sjr{s)) 
11 



(27) 



n{\k 


pP ') 




p 9a) 



where 

fcg/\{0} 

Applying the inverse Laplace transformation to (l27|), we find S{t). 

The numerical results follows. The typical curves S{t) within a certain 

time window are shown in figure 5 and 6. It is clearly seen that there are three 
characteristic kinetic modes. At high temperatures (curve 1 in Fig 5), the kinetics 
S{t) is governed by the exponential law. This mode is observed in the case of 
Teq ^ where T^q is the characteristic time during which the distribution f{x,t) 
reaches quasi-equilibrium, and A is the reaction sink parameter. In this case, during 
the observation periods t ^ r^g, the function f{x,t) is close to a homogeneous 
distribution and the population of the reaction sink area (disc Zp) is equal to p~^S{t). 
Therefore, S{t) is determined by the kinetic equation 

which describes the exponential relaxation S{t) = exp(— t/r) with the relaxation 
time r = p^ /X. Since the limiting kinetic stage here is the reaction sink, this mode 
may be called the reaction control mode. As the temperature decreases (curves 2 
and 3 in figure 5), a small section of power relaxation appears on the curve S{t) and 
then extends to the major part of the time window of observation. This change in 
the curves S{t) corresponds to the transition regime for which, with the decrease of 
temperature, the time of relaxation to equilibrium, Tgg, and the time parameter of 
the reaction sink become commensurable. In this case, the kinetics of S{t) is 
limited first by the ultrametric diffusion (power section of the curve S{t)), and then 
by the reaction sink (the exponential decay). Note that if the time window under 
observation shows only the power section of the curves S(t), it seems that the decay 
rate S(t) is growing with the decrease of temperature. 

With further decrease of temperature (curves 1-3 in figure 6), a beak-shaped 
section appears on the curve S{t) which then extends to the entire time window 
under observation. This mode corresponds to situations with r^g ^ In this case, 
the decay rate S{t) is limited by filling of the reaction sink area due to ultrametric 
diffusion. It is natural to call such a mode the diffusion control mode. Now, there 
are no "anomalous" kinetic effects in the time window under observation, the rate 
of decay S{t) falls down as the temperature decreases. 

It is interesting to note that all peculiarities of the kinetic curves for the model 
considered above are also typical for the ligand rebinding kinetics of myoglobin 



0, |2^. Naturally, this gives grounds for some optimism. Nevertheless, we would 
like to stress that the main objective of this paper is to demonstrate a method for the 
construction of ultrametric diffusion models based on p-adic analysis. The examples 
given above should be considered as a mere illustration of the possibilities offered 
by this approach, rather than particular models giving a quantitative description of 
particular experiments. The development of such models seems to be a special task 
in each particular case. 
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time (in dimensionless units) 



initial distribution 
diffusion ■ I diffusion 
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time (in dimensionless units) 
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